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An analytical model for calculation of ascent steady-state tile loading was 
developed and validated with wind-tunnel data. The analytical model is described 
and results are given. Results are given for loading due to shocks and skin fric- 
tion. The analysis included calculation of internal flaw (porous media flow and 
channel flow) to obtain pressures and integration of the pressures to obtain forces 
and moments on an insulation tile. A heat- transfer program was modified by using 
analogies between heat transfer and fluid flow so that it could be used for internal- 
flow calculation. The type of insulation tile considered was undensified reusable 
surface insulation (RSI) without gap fillers, and the location studied was the lower 
surface of the orbiter. Force and moment results are reported for parameter varia- 
tions on surface pressure distribution, gap sizes, insulation permeability, and tile 
thickness. 


INTRODUCTION 

The NASA Structures Team was established in February 1980 for the purpose of 
resolving problems of the Shuttle Thermal Protection System (TPS) as they existed at 
that time. This work was performed to support the NASA Structures Team in the inde- 
pendent structural assessment of the TPS for the first launch of the Space Transpor- 
tation System (STS-1). The present analysis of TPS steady aerodynamic ascent loads 
was needed because the existing analysis did not agree with wind-tunnel data. The 
purpose of this report is to describe the detailed analysis which was done at the 
Langley Research Center (LaRC) and to present some selected results. The analysis 
was accomplished by formulating a mathematical model for calculations of steady aero- 
dynamic loads by adapting the Martin Marietta Interactive Thermal Analysis System to 
perform internal steady flow analysis. The M1TAS (ref. 1) program could be used for 
internal- flow calculation by modifying the input to the program based on analogies 
between heat transfer and fluid flow. The program was used to calculate internal 
mass flows and pressures of an array of nine LI900 undensified tiles (9 lb/ft^ 
density). A postprocessor loads calculation program was written to calculate loads 
using output from the MITAS program and other inputs. This model was used to 
calculate tile loads caused by stationary shocks. 

This report consists of sections on flow modeling and validation, force and 
moment postprocessor, and results of parametric runs. A program user guide and a 
copy of the MITAS input computer code of the model are available in reference 2. 

This computer model was subsequently modified for use in gap heating analysis 


(described 

in ref. 3). 
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O 

area on which pressure difference acts, in 

p 

area which is stressed, in 
discharge coefficient 
local skin- friction coefficient 
tile thickness, in. 
force, Ibf 

force in the x-direction, lbf 
force in the y-direction, lbf 
force in the z-direction, lbf 

conductance for heat transfer, Btu/hr-°F; conductance for fluid flow, 
lbm/sec-psi 

2 

mass unit conversion, lbm-f t/lbf-sec 

length of moment arm, in. 

permeability, i n“ 

tile roughness factor, in. 

spring constant , lbf/in. 

thermal conductivity, Btu/hr-f t-°F 

length of tile diagonal, 8.49 in. 

SIP height after tile installation, in. 

filler-bar thickness after tile installation, in. 

moment, in- lbf 

mass-flow rate, lbm/sec 

pressure, psi 

pressure difference, psi 

pressure difference between surface and SIP, psi 
pressure difference between tile gap and filler bar, psi 
pressure at bottom of gap, psi 
pressure at top of gap, psi 
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pressure on tile surface, psi 
pressure at tile bond line, psi 
heat-flow rate, Btu/hr 

distance from center of one node to next, in, 

temperature, °R 

tile bond line 

free- stream velocity , fps 

tile element width, in. 

width of gap between filler bar and tile 

minimum gap between filler bar and tile 

initial gap width 

width of gap between tiles 

tile width, in. 

Shuttle axis system 
local (tile)axis system 

distance from airfoil leading edge, in. (fig. 13) 
stream wise integration step, in. 
distance from shock , x ~ x ^ in. (fig. 13) 
local (tile) coordinate system 

distance from airfoil leading edge to tile leading edge, in. (fig. 13) 

distance from tile leading edge to foot of shock, in. 

distance from airfoil leading edge to shock location, in. (fig. 13) 

distance along Y-axis with respect to Shuttle coordinate system 

boundary-layer thickness, in. 

viscosity, lbm/in-sec 

viscosity of free stream, lbm/ft-sec 
density, lbm/in3 

free-stream density, lbm/ft^ 
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2 

T w surface shear stress unaccounted for shock, lbf/ft 

2 

T surface shear stress accounted for shock, lbf/ft 

Subscripts: 


i/j tile matrix location 

\ local (tile) 

X x-direction 

Y y-direction 

Z z-direction 


INTERNAL FLOW MODEL USING MTTAS 

A three-dimensional analytical model of an array of nine LI900 undensified tiles 
was developed to adequately represent the internal flow through the Shuttle thermal 
protection system. The primary requirement was to obtain pressures in the tile and 
in the gaps between tiles in sufficient detail to allow tile loading to be deter- 
mined. MITAS (ref. 1) was used in preference to writing a custom computer program to 
save time and manpower. MITAS (Martin Marietta Interactive Thermal Analysis System) 
is a general purpose finite-difference heat- transfer program and was applicable to 
this analysis of internal incompressible (low subsonic Mach number) flow because of 
the similarity in the equations for internal incompressible flow and the equations 
for heat transfer. 

The first step in the modeling was to establish in detail how the MITAS program 
would be used to model incompressible internal flow. The Laplace equation is appli- 
cable to both incompressible flow and heat transfer. These types of problems can be 
converted to an analog representation consisting of nodes connected by conductors, 
which is precisely the way in which MITAS is used to represent a heat- transfer prob- 
lem. This type of representation “ commonly called a lumped node, electrical analog 
representation of a heat- transfer problem “ also applies to internal incompressible 
flow. The required input to the MITAS program is the description of nodes and con- 
ductors for the electrical analog representation of the thermal or fluid flow prob- 
lem. MTTAS could be used directly for this flow problem because of the exact 
correspondence between pressure and temperature and between mass flow and heat 
flaw, Some comparisons of conductors for heat flow and fluid flow are given in the 
appendix. 

The input to the MITAS program for the steady-state flow situation is first of 
all the node and conductor definition for the electrical analog representation of the 
geometric configuration. The network and corresponding MITAS input was first devel- 
oped for a single tile; then with the single-tile representation as a guide, the 
network for the multiple - tile model was developed. Three-dimensional modeling was 
required to obtain internal pressure for tile, strain isolation pad (SIP), and 
filler- bar and tile- to- tile gap pressures from given surface pressure boundary 
conditions. 

A single- tile three-dimensional model was developed first. A sketch showing the 
flow paths through a single tile is shown in figure 1» The ceramic coating which 
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covers most of the tile surface is assumed t o be impermeable. Flow through cracks 
and chipped areas was neglected. The flow through the silicone rubber bond line on 
the bottom of the tile was also considered t o be insignificant. Thus, areas where 
flow can enter or leave a tile are the terminator gap area and the 1/2-in-wide 
unbonded area on the bottom outer edge of the tile. Flow can enter or leave the SIP 
through the unbonded area on the bottom of the tile or the gap between the tile bot- 
tom and the filler bar. The connections to the surface through which all internal 
flow enters or leaves are the gaps between the tiles. With the flow paths in the TPS 
established, the next step was to define the nodes and conductors for the single- 
tile model. The dimensions of the single- tile model are 6 by 6 in., and the thick- 
ness of the tile and SIP are variable. The SIP thickness was set at 0.16 in. in 
this analysis. The nodal arrangement for the subject tile is shown in figures 2 
through 4. Nodes 61 through 96, which are the nodes inside the tile at the surface, 
are shown in figure 2. Nodes 25 through 60 (inside the tile at the bond line plane), 
nodes GI through G24 (at the bottom of the gap between tiles), and nodes B201 through 
B224 (the boundary nodes for the internal flow which are located at the surface over 
the gap between tiles) are shewn in figure 3. Nodes 97 through 132, which are in the 

SIP, are shown in figure 4. Conductors were defined to connect the nodes in three 

dimensions. Examples of nodes in the x^-y ^ plane that were connected by conductors 
would be 61 to 62, 61 t o 67, 25 t o 26, 25 to 31, 97 t o 98, and 97 to 103. Examples 
of nodes that were connected by conductors for flow in the z^^-direction would be 61 
to 25, 62 t o 26, 67 t o 31, 68 t o 32, 25 t o 97, 26 t o 98, 31 to 103, and 32 t o 104. 

There are two parallel flow paths from the bottom of the gap to the bottom of the 

tile: one flow path through the edge of the tile and one flow path under the edge of 

the tile between the tile and filler bar. (See fig. 1, view A.) These flow paths 
are represented in figure 3 with electrical resistance symbols connecting gap nodes 
to tile bond line nodes such as nodes G2 and 26, for example. Figure 3 also shows 
the conductors connecting surface boundary nodes to nodes at the bottom of the gap, 
such as the conductor shown connecting boundary node B202 to gap node G2. 

With the single- tile model as a guide, a model of an array of nine tiles was 
developed by adding eight more similar tiles to the single- tile model. Flow from one 
tile to another is through the filler bar and through the gaps between tiles. The 
overall pattern of tiles for the multitile model is shown in figure 5. A staggered 
pattern was chosen because this is the dominant pattern of tiles on the Shuttle. The 
nine tiles in the model are shown numbered 1 through 9. Tile 9 is the subject tile, 
also referred to as the "0000" tile in figures 6, 7, and 8. The locations of the 
nodes in the multiple- tile array are shown in figures 6, 7, and 8. Figure 6 shows 
nodes inside the tile at the surface for each of the nine tiles. Figure 7 shows 
nodes at the tile bond line and at the surface and bottom of the adjacent gaps for 
each of the tiles. Figure 8 shows nodes in the SIP for each of the tiles. 


Modeling of Conductance 

After defining the nodes and conductor connections for the nine- tile model, the 
conductance values for these conductors were determined. The general categories of 
conductance values t o be determined were tile material, SIP material, gap between 
tiles, tile terminator gap, gap between tile and filler bar, and gap between filler 
bar and SIP. 
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The conductance for flow through tile and SIP material was based on flow 
through porous media (ref. 4). The conductance for flow through porous media 
as used in the analysis is determined as follows: 


G 


K pA 

P c 
p As 


The flow through the tile in the terminator gap location, the tile bond line 
node at the edge of the tile (node 55 i n fig. 3, for example), was modeled as flow 
through porous media by using an equivalent flow area. This equivalent flow area was 
the average of the terminator gap area (the product of terminator gap width and node 
width) and the node cross-sectional area for flow (the product of D . /2 and node 

width). tllS 


The equations which represent laminar flow between parallel walls were used to 
represent flow in the gap between tiles. The validity of using these equations was 
demonstrated in reference 5 by comparison of the theoretical results for laminar flow 
between parallel walls to experimental data for flow in the gaps between tiles. The 
equations for mass flow and conductance, from reference 6, for laminar flow between 
parallel walls as used in the analysis are as follows: 


m = - 


. PC A n 

1 t c dp 


1 2 


p ds 


, pwi> 

1 r t c 
1 2 p As 


tile - to - tile gap width 
flow area, W^. x Gap length 

pressure gradient 

distance from center of one node to next 

As can be seen from the equation for flow conductance, the conductance varies as the 
cube of the gap width. Consequently, it was necessary to account for changes in gap 
width caused by forces on the tiles for the gap between tiles and the gap between 
filler bar and tile. 

Tire gap size, and resulting conductance for flow through the gap between tiles, 
was calculated iteratively for tiles which had significant size forces on them due to 
gap pressures. These side forces result from the surface pressure gradient, such as 
the gradient due to a shock being positioned over a tile. The tiles which required 
gap variations to be considered were the tiles which had the shock positioned on them 
coinciding with the shock being positioned on the subject tile. Figure 9 shows the 
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gaps between tiles in the global model, which were allowed to vary. The calculations 
of tile deflection were done only on the subject tile, tile 9, and the gap deflec- 
tions of the other tiles, tile 3 and tile 7, were set to the same values as the cor- 
responding locations on the subject tile. All other gaps between tiles were fixed at 
the initial gap opening. Gap pressures on the sides of the subject tile were inte- 
grated to obtain a resultant side force in the x - and y^-directions. Rotation of 
the tile about the 2, -axis being neglected, the translational motion of the subject 
tile was then obtained from a spring representation of the SIP according to the fol- 
lowing equations : 


Ax = 


X 

K 


Ay, = 


K 


where 

Ax 

t 

F X 

Ay, 

f y 

K 


subject tile x - deflection 


force in the x,-direction 


subject tile y -deflection 
force in the y -direction 
lateral spring constant for SIP 


The SIP lateral spring constant used in this analysis was 110 lbf/in» (ref. 7). The 
gaps shown solid in figure 9 decrease from the initial value by the tile deflection, 
whereas the gaps shown hatched increase. A limiting condition for gap closure is the 
effective gap which exists when the two nonsmooth surfaces are in contact. The gaps 
between tiles were adjusted based on the calculated pressures until the calculated 
gap sizes converged. 


The gap size and resulting conductance for flow through the gap between tile and 
filler bar was calculated iteratively for all tiles in the global model. The size of 
this gap is determined by the relative normal deflection between the bottom surface 
of the tile edge and filler bar and also by a silicone rubber membrane effect. The 
size of this gap is determined by the pressure forces which act on the filler- bar 
membrane and tile, causing deflections from the initial position. 


The dimensions and deflections which influence this gap are shewn in figure 10. 
The equation which was used for this tile edge deflection due to stretching or com- 
pressing the SIP is 


Ap,A 

AJL = ?-?- ■ JL 
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where 


Ap^ pressure difference between external surface and SIP 

A i area on which Ap^ acts, loading SIP 

A^-j area of SIP which is stressed 

SIP thickness after tile installation ( 0.1 6 > in. ) 

E 1 modulus of elasticity for SIP 

The filler-bar deflection is caused by the pressure differential across the silicone 
rubber surface membrane and is defined by 


a P 9 A 

hi 2_p2_ ^ 

2 E A ^ X 2 

2 02 


where 

Ap 2 pressure difference between tile gap and filler bar 

Ap2 area of filler bar on which Ap2 acts 

A ft1 area of filler bar which is stressed 

^2 filler-bar thickness after tile installation 

modulus of elasticity for filler bar 

The desired width of the gap between the filler bar (FB) and tile is obtained based 
on the following equation: 

= AJL “A A, + W 
r z t o 

where 

Wf width of gap between filler bar and tile 

W Q initial gap width 

A limiting condition for this gap width is 


W 


f 


> w 


min 


where W . is the effective gap which exists when the tile edge touches the filler 
nun oi- o 

bar (0.001 in.). The final condition which was considered l n the determination of 
the gap size was the silicone rubber membrane effect. This effect causes a check 
valve action which restricts flow out of the SIP/FB more than it restricts inflow. 
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It was assumed that this membrane effect would dominate the gap size determination 
such that for this analysis it was adequate to assume that the gap height for flow 
out of the SIP/PB would be and the gap height for flow into the SIP/FB would 

be as calculated from the aboflffe n equation for w The gaps between the tile and 
filler bar and the resulting flow conductance bSsed on laminar flow between parallel 
walls were also changed based on the calculated pressures until the calculated gaps 
converged. 

Flow through the gap between the SIP and filler bar was modeled as flow through 
a slot. The discharge coefficient for this slot was assumed to be 0.5. The 

flow crpss- sectional area for tiles with significant side deflections was varied 
consistent with the side deflections calculated for the deflections of the gap 
between tiles mentioned previously. 


Discussion of Parameters Which Influence Flow i n TPS 

The parameters which influence flow in TPS can be grouped into the following 
categories : 

1. Surface pressure distribution 

2. Tile- to- tile gaps 

/ 

3. Terminator gaps 

4 . Internal tile flow 

5. Tile-to-SIP flow path 

6. Gap between tile and filler bar 

7. Internal SIP flow 

8. Gap between SIP and filler bar 

In each of these categories, all known parameters were evaluated to determine which 
ones were significant enough t o be included in the flow model. 

The first area to be discussed is the surface pressure distribution. The param- 
eters which were considered for this category in the model in order of importance 
were shock strength, shock location, shock smearing, absolute pressure, and absolute 
temperature. Shock location was defined as the location of the foot of the shock on 
the tile, and shock strength was defined as the pressure difference between the foot 
and the peak. Shock smearing has t o do with the distance over which the shock is 
spread, which is normally governed by boundary-layer thickness. Absolute temperature 
influences viscosity and density of gas, whereas absolute pressure produces mainly a 
density effect. 

In the category of tile-to-tile gaps, the parameters considered were minimum gap 
width, initial gap width, SIP shear" stiffness, tile thickness, and gap pattern. The 
initial gap occurs in the unloaded condition, whereas the minimum gap is the gap 
which represents tile - to - tile contact. SIP shear stiffness affects the side-force 
deflection characteristics of the tile; these side deflections are important in 
determining gap variations. The tile thickness and tile pattern affect the geometry 
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of these flow paths. The model was limited to one gap pattern: a staggered pattern 

of square tiles. 

In the category of the terminator gap, the only parameter considered was the 
terminator gap width. 

In the category of internal tile flow, the parameters considered were tile 
thickness and tile permeability ■ 

In the category of tile-to-SIP flow path, the parameters considered were effect 
of tile side deflection on the flow path, dimensional tolerancing on the SIP and 
filler bar, and the tile to SIP discharge coefficient. The dimensional tolerancing, 
side deflections, and discharge coefficient describe the flow through the slot 
between the SIP and filler bar. 

In the category of flow through the gap between the tile and filler bar, the 
parameters considered were membrane effects of the filler- bar silicone coating, ini- 
tial gap size, SIP and filler- bar stiffness, and minimum gap effect. The membrane 
effect is a tendancy of the edge of the filler-bar coating to roll up against the 
bottom of the tile producing a checking effect on flow in the direction from the SIP 
to the gap between tiles. The same initial gap size and minimum gap size were spec- 
ified for all tiles in the global model. The minimum gap size represents the con- 
dition where the edge of the tile is in contact with the filler- bar surface membrane. 
The stiffnesses of the SIP and filler bar were represented with an effective modulus 
of elasticity. 

In the category of internal SIP flow, the parameters considered were SIP 
permeability, and SIP deflections. The permeability of the SIP was obtained from 
reference 6. These data had t o be adjusted to account for differences between test 
conditions and expected flight conditions (the SIP expands under load in flight). 

In the category of the gap between SIP and filler bar, which is a gap between 
the silicone membranes of the SIP and filler bar, the parameters considered were a 
discharge coefficient, the effect of tile side deflection, and the effect of dimen- 
sional tolerance. 

The more significant categories in the model were surface pressure distribution, 
the gap between tile and filler bar, tile-to-tile gaps, and internal SIP flow. Less 
significant were the terminator gaps, internal tile flow, and tile-to-SIP flow path. 
The least significant category considered in the model was the gap between SIP and 
filler bar. 

There are many other parameters influencing flow in TPS which were not con- 
sidered significant enough to include in the flow model. Some relevant parameters 
which were not considered in the flow model, listed by category, are as follows: 

1. Surface pressure distribution: shock stability and tile surface step effect 

2. Tile-to-tile gaps: compressibility effects in the gaps, influence of tile 

rotations on gap size, and initial gap variation from tile to tile 

3. Terminator gap: changes in tile porosity due to the coating application and 

gap machining process and variation of gaps from tile to tile 
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4. Internal tile flow: permeability of tile coating and bond line due to 

imperfections 

5. Tile-to-SIP flow path: the effect of rotational deflections and 

z ^-deflections of the tile 

6. Gap between tile and filler bar: condition of the filler- bar transfer coat 

edge, condition of the tile bottom surface, changes in filler bar and SIP 
caused by proof testing, and dimensional tolerance from tile to tile 

7. Internal SIP flow: effect of proof testing, directional bias of SIP, and 

wicking of silicone rubber into the SIP 

8. Gap between SIP and filler bar: dimensional variations from tile to tile and 

changes due to the bonding process 


FORCE AND MOMENT MODEL 

The force and moment model is a postprocessor routine of the overall MITAS tile 
flow analysis model. The force and moment model numerically integrates the pressure 
profiles produced by the MITAS model to determine the forces and corresponding 
moments on the Shuttle tiles. The model accounts for the pressure distributions 
across the tile surface, at the tile bond line, and at the top and bottom of the gap 
around the tile perimeter. The model also includes skin- friction effects . The 
analysis employs a right-hand coordinate system and resolves the forces and moments 
in the x-, y-, and z-directions (fig. 11). 

The tile was separated into 1-in elements producing a 6 x 6 matrix on the tile 
surface and at the tile bond line, and 24 elements in the gap around the tile perim- 
eter (fig. 12). The pressures on the tile surface and at the tile bond line were 
assumed t o be constant per element. In the gap, the pressure was assumed to follow a 
linear distribution from the top to the bottom of the gap at each station around the 
tile. The pressures on the tile surface and at the tile bond line contribute to the 
normal forces and the moments about the X- and Y-axes. The gap pressures around the 
tile perimeter contribute to the forces in the x- and y-directions and the moments 
about the X-, Y-, and Z-axes. The skin friction contributes to the force in the 

x-direction and moment about the Y-axis. The calculated forces and moments are added 
together to obtain the resultant vectors, F , Fy, and F^ and the resultant 
moments about the X-, Y-, and Z-axes, My, and Mg. 

Normal Forces and Moments 

With the use of the pressure distributions on the tile surface and tile bond 
line, the normal force „ acting on the tile is expressed per element by 

*Z 


F Z. . ( P TBL. . P sur..\ A el.. 
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Therefore, the resultant normal 


where 1 < i < 6, 1 < j < 6, and A 
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The moments about the X- and Y-axes produced by the normal forces and summed per 
element are defined by 
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where y . . = (i - j + 6) (0.7071) and x • . 


(i + j - 1) (0.7071) . 


Gap Forces and Moments 

In the gap around the tile perimeter, the perpendicular force on each elemental 
face, the moment, and the moment arm are defined by 
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where 1 < i < 24. The components in the x- and y- direct ^ ons are defined by 
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where 1 < i < 24 and 


S„ = +1 


S,, = +1 


S = -1 


(1 < i < 6 and 19 < i < 24) 

(7 < i < 18) 

(13 < i < 24) 

(1 < i < 12) 


The moment about the Z-axis per element is expressed by 



F i L i 


where 
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and 1 < i < 24. The resultant moment about the Z-axis due to the gap pressures is 
determined by the sum 


M 


Z 




Skin- Friction Forces and Moments 

The skin- friction contribution to the force and moment calculation is based on 
the aerodynamic loading on the tile. The local shear stress x^ on the tile 
surface is determined by the typical equation 

x = CLp U 2 
w 2a r 00 00 
c 


where is the local skin- friction coefficient. The skin- friction coefficient is 

assumed to follow the resistance formula for a uniformly rough flat plate (ref. 6): 


_ (2,87 + 1 . 58 log x 

f = l 

2 6 

which is valid when 10 < x/K r < 10 , (See fig. 13.) The shear stress is then 

adjusted to ~ to account for local boundary- layer separation directly beneath a 
shock on a til^f. This adjustment is accomplished by determining a shear- stress 
reduction factor (x x ) as a function of a shock location and the boundary-layer 
thickness (ref. 8). The boundary-layer thickness 6 is estimated by the smooth 
plate equation (ref. 6): 


6 = 


-0.37x1 


^oU„x/12j 


-1 .5 


Given the free- stream conditions (i.e., 

and the running length x where x < x < x + L, the 

so that the factor ^/x can be found (ref? 8) . The 

calculated and adjusted to x w to account for a shock 

area on which x„ acts is determined by 
w J 


and U ), the shock location 

CO 

ratio x'/S can be determined 
shear stress, x w is then 
on a tile. The differential 
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The corresponding differential forces and moments are then determined by 
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% AA 
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144 


AM y - AF X D tile 


Thus, the forces and moments due to the skin friction are 


F = y AF 
X ^ X 

M y = E ^ 


RESULTS AND DISCUSSION 

The results of the analysis which are presented herein are a validation of the 
model and the calculated loads for a range of parameters. Tbs model was validated by 
comparing calculated pressures with measured pressures from wind-tunnel tests. The 
load sensitivities of eight different parameters were determined by calculating 
forces and moments for varying input data to the model. 


Comparison With Test Data 

The pressure calculations of the flow model were validated by comparison with 
test data. The OS36 (fig. 14) and OS55 (fig. 15) Shuttle tile wind-tunnel test data 
were used for comparison with calculations (refs. 10 and 11, respectively). Ihe "OS" 
in the test designation stands for orbiter structural test. For the OS36 and OS55 
tests, arrays of tiles mounted on SIP and equipped with pressure measurement instru- 
mentation were placed in a wind tunnel and subjected to aerodynamic shock environ- 
ments. The purpose of these tests was to observe damage and measure internal 
pressures caused by shocks. Figure 16 shows comparison between calculated and mea- 
sured pressures at the tile bond line at the diagonal centerline for the OS36 wind- 
tunnel test. The pressure at the tile bond line is the pressure inside the tile 
adjacent to the silicone rubber bond which attaches the tile to SIP. The rectangle 
at the bottom of the plot indicates the location of the 6- by 6-in. tile which is at 
a 45° angle so that the overall length in the flow direction is 8.49 in. Comparison 
of calculated and measured pressures at the tile bond line for the OS55 wind-tunnel 
test, runs 4:4, 4:5, and 4:6 are presented in figures 17 , 18, and 19. The results of 
the calculations were in good agreement with the measured data; therefore, the model 
could be used to calculate forces and moments on tiles for STS-1 and no additional 
wind-tunnel testing of tiles was needed. 
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Parametric Runs of Multi ti le Model 


A series of runs of the model was made to determine the influence of tile param- 
eters on tile steady aerodynamic loads due to shocks. The tile parameters which were 
examined for sensitivity with the aid of the model are 


1. Initial tile - to - tile gap 

2 Terminator gap width 

3. Gap between filler bar and tile 

4. SIP permeability 

5. Tile permeability 

6. Tile thickness 

7. Pressure effect (altitude effect) 

8. Boundary-layer thickness 

Aerodynamic force and moment calculations were made for several different shock 
pressure profiles over a range of shock locations from the leading to the trailing 
edge of the subject tile. 

The parameters were studied with estimated pressure profiles for the mid- 
fuselage zone 6 (MF-6) region of the Shuttle orbiter (fig. 20) and measured pressure 
profiles from the OS55 wind-tunnel test. The pressure profiles and results of the 
parametric analysis are shown in figures 21 through 27. Figure 21 shows estimated 
MF-6 pressure profiles for shocks of 2 psi and boundary-layer thicknesses of 2.5 and 
5.0 in. High and low upstream pressures were considered, high being 5.4 psi and low 
being 1.8 psi. Figure 22 shows the forces and moments produced by the pressure 
profiles of figure 21 for tile thicknesses of 1.0 and 2.0 in. and for various initial 
tile-to-tile gaps. Estimated MF-6 pressure profiles for two specific tiles in this 
region for the STS-1 ascent condition are shown in figure 23. The first tile is tile 
394020-272, 1.47 in. thick, and is referred to as tile A. The second tile is tile 
394034-087, 1.0 in. thick, and is referred to as tile B. Forces and moments for a 

nominal configuration of tile A and tile B are shown in figure 24, and the effect of 
initial gap variations on force and moment is shown in figure 25. The calculated 
forces and moments for tile A for single parameter variations of terminator gap, gap 
between filler bar and tile, SIP permeability, and tile permeability are presented in 
figure 26. The calculated forces and moments for tile B for single parameter varia- 
tions of terminator gap and filler bar to tile gap are shown in figure 27. The 
measured pressure profiles from the OS55 wind-tunnel test are shown in figure 28. 

The forces and moments produced by the OS55 pressure profile for single parameter 
variations of initial tile- to- tile gap and tile thickness are shown in figures 29 
and 30. 

Several of the parameters studied were found to have a very significant effect 
on the tile loading. The location of the foot of the shock and the steepness of the 
shock were found to be significant in all cases where these parameters were varied. 

An increase in upstream pressure from 1.8 to 5.4 psi, an increase of 3.6 psi in the 
pressure distribution, caused the maximum forces and moments to increase by approxi- 



mately 50 percent as shown in figures 22(a) and (b). When the initial tile - to - tile 

gap was decreased from the nominal value of 0.050 to 0,015 in., the maximum normal 

force for the MF-6 high pressure profile was decreased by 30 percent as shown in 

figure 22(c). The separated shock cases plotted in figure 25 and OS55 pressure 

profile results plotted in figure 29 also shew significant effects when the tile-to- 

tile gap is decreased from the nominal value of 0.05 in. SIP permeability has a 

strong effect on normal force, as shown in figure 26(e), when the SIP permeability 

was increased by a factor of 4 for the separated shock, MF-6 profile, the maximum 

normal force increased by a factor of 2. Tile thickness has a significant effect on 

loading as shown in figure 30. The highest combined loading is calculated for the 

case where the foot of the shock is located at x /L = 0.73. 

s' 

CONCLUDING REMARKS 

A model was developed to calculate steady aerodynamic loads for Shuttle tiles. 
The use of the M1TAS heat- transfer program to calculate internal incompressible flow, 
flow through porous media, and channel flow was demonstrated. Calculated pressures 
at the tile bond line caused by aerodynamic shock on the tile surface were in good 
agreement with measured pressures at the tile bond line from wind-tunnel tests. The 
model was used to calculate forces and moments on TPS tiles for STS-1. 

The results of the analysis were used by the NASA Structures Team to obtain tile 
loadings for the STS-1 flight of the Shuttle for conditions where wind-tunnel data 
were not available. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 26, 1984 
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APPENDIX 


EQUATIONS FOR HEAT TRANSFER AND FLUID FLOW 


For heat flow. 


q = kA 


dT 
c ds 


kA 


G = 


As 


where 

q 

k 

dT 

ds 


heat flux 

thermal conductivity 
temperature gradient 


A cross-sectional area for flow 

c 

As distance from center of one node to next 

For flow through porous media (ref. 11), 


m = 


K pA 

c dP 
p ds 


G 


K pA 
_E £ 

P As 


where 

m 


p 

p 

dP 

ds 


mass-flow rate 
permeability 
viscosity 
density 

pressure gradient 
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APPENDIX 


For flow through slots and orifices, 

m = CdAc(2pg /Ap) 1 ^ 2 Ap 
c 


G = C A ( 2pg /Ap) 
a c c 


1/2 


where 

m mass-flow rate 

discharge coefficient 
A c flow area 

P density 

g mass unit conversion 

pressure difference 

Units which are commonly used for heat-flow rate and for temperature can also be 
compared with the units for the corresponding mass-flow rate and pressure which were 
chosen for the fluid flow model. For the thermal model, common units are 

• 

<3 heat-flow rate, Btu/hr 

T temperature, °P 

G conductance, Btu/hr-°P 

For the fluid flow model, the units chosen were 
m mass-flow rate, Ibm/sec 

P pressure, psi 


G 


conductance, Ibm/sec-psi 
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Figure 6— Tile array showing surface nodes. 
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Figure 7.- Tile array showing nodes at bond line. 
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Figure 9.— Tile array showing gap changes due to x - and y -deflections 
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Figure 21.- MF—6 pressure profiles. 
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Figure 23.— MF-6 pressure profiles . 
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Figure 27.- Calculated normal force and moment as function of shock location for MF-6 region, pressure 
profile B, separated shock, LI900 undensified tile 394034-087, and tile thickness of 1 in. 



Terminator gap, in 
0.045 



Figure 27.- Continued 
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Figure 30.- Calculated normal force and moment as flection of thicker ss for 

OS55 pressure profile. 




Thickness, in. 


4 





1. Report No. 

NASA TP-2255 


2. Government Accession No. 


3. Recipient's Catalog No. 


4. Title and Subtitle 

5. Report Date 

STEADY INTERNAL FLOW AND AERODYNAMIC LOADS 

March 1984 

ANALYSIS OF SHUTTLE THERMAL PROTECTION 

6. Performing Organization Code 

SYSTEM 

505-31-53-08 

7. Author(s) 

8. Performing Organization Report No. 

Dennis H. Petley, William Alexander, Jr., 
George W. Ivey, Jr., and Patricia A. Kerr 

L- 15695 

10. Work Unit No. 

9. Performing Organization Name and Address 


NASA Langley Research Center 
Hampton, VA 23665 

11. Contract or Grant No. 


13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

Technical Paper 


14. Sponsoring Agency Code 

15. Supplementary Notes 


16. Abstract 


An analytical model for calculation of ascent steady-state tile loading was 
developed and validated with wind-tunnel data. The analytical model is described 
and results are given. Results are given for loading due to shocks and skin fric- 
tion. The analysis included calculation of internal flow (porous media flow and 
channel flow) to obtain pressures and integration of the pressures to obtain forces 
and moments on an insulation tile. A heat- transfer program was modified by using 
analogies between heat transfer and fluid flow so that it could be used for internal- 
flow calculation. The type of insulation tile considered was undensified reusable 
surface insulation (RSI) without gap fillers, and the location studied was the lower 
surface of the orbiter. Force and moment results are reported for parameter varia- 
tions on surface pressure distribution, gap sizes, insulation permeability, and tile 
thickness. 


17. Key Words (Suggested by Author(s)) 

Thermal Protection System 

Ascent loading from shock flow 
Internal incompressible flow 
Porous media 
Channel flow 


18. Distribution Statement 

Unclassified 


19. Security Classif. (of this report) 

Unclassified 


20. Security Classif. (of this page) 21 . No. of 

Unclassified 67 


" Unlimited 


Subject Category 18 

22. Price 

AO 4 


For sale Jby the National Technical information Service, Springfield, Virginia 22161 


NASA-Langley, 1984 

















